Origin of large thermopower in LiRh 2 04 
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Motivated by the newly synthesized mixed-valent spinel LiRl^CU for which a large thermopower 
is observed in the metallic cubic phase above 230K [Okamoto et al, Phys. Rev. Lett. 101, 
086404(2008)] , we calculate the Seebeck coefficient by the combination of local density approximation 
and dynamical mean field theory (LDA+DMFT). The experimental values are well reproduced not 
only by LDA+DMFT but also by the less involved Boltzmann equation approach. A careful analysis 
of the latter shows unexpectedly that the origin of the large thermopower shares a common root with 
a very different oxide: Na x Co02. We also discuss how it is possible to further increase the power 
factor of LiRti2 04 through doping, which makes the material even more promoising for technological 
applications. 



I. INTRODUCTION 

Designing and searching for good thermoelectric ma- 
terials have a long history of extensive studies due to the 
scientific interest and potential technological importance, 
particularly for generating electrical power from heat 
(gradients) and for cooling through the Peltier effect^. 
Hitherto, the main target materials have been various in- 
sulators or semiconductors such as B'^Te^ and FeSb2 2 , 
since it was believed that huge thermopowers cannot be 
expected for metals. However, recently, novel metallic 
systems with large thermopower have been discovered 
and attracted much attention. Generally, materials with 
strong electronic correlations are promising 3 ; and a fa- 
mous example is Na^CoO^, for which a metallic resistiv- 
ity as low as p=0.2 mf2cm and a thermopower as large 
as ,5=100 £tV/K are observed simultaneously at 300K 4 . 
The coexistence of low resistivity and large thermopower 
results in a large power factor (S 2 /p), which is especially 
important for device applications. 

Most recently, Okamoto et al. synthesized a new 
mixed-valent spinel oxide, LiRli204 5 . This novel oxide 
shows two structural phase transitions, i.e., the cubic- 
to-tetragonal transition at 230K and the tetragonal-to- 
orthorhombic transition at 170K. Particularly interesting 
is however the high temperature cubic phase: Despite the 
metalicity which is reflected in a small resistivity and the 
existence of a Fermi-edge, the thermopower is as large as 
80/iV/K at 800K, which is exceptional for metallic sys- 
tems. 

On the theoretical side, a variety of studies have been 
performed to understand the mechanism of large ther- 
mopowers in metallic systems. Among others, Koshibae 
et al. derived an expression for the Seebeck coefficient 
of strongly correlated systems in the high-temperature 
limit 6 . Considering the orbital and spin degrees of free- 
dom of localized electrons, they estimated the thermo- 
electric power of Naa;Co02 to be 150 fj,V/K. 

However, when the temperature (T) is much lower than 



the energy scale of the band width (~ 2 eV), it is ex- 
pected that the band dispersion of the system also plays 
a crucial role as has been suggested from first princi- 
ple (bandstructure) studies^. Indeed, recently, two of 
the present authors proposed that the peculiar shape of 
the valence band (the so-called a\ g band) is important 
to realize a large thermopower and high conductivity in 
Na x Co02 9 . The different theoretical proposals led to a 
heated discussion^ and also to a proposal to discrimi- 
nate between them through the respective temperature 
dependence^. 

The motivation of the present study is to clarify the 
origin of the large thermopower in LiRli2 04. For this 
purpose, we first perform a LDA+DMFT— calculation 
(the combination of the local density approximation and 
the dynamical mean field theoryi 3 -) , employing the Kubo 

formula for the Seebeck coefficient^. This ah initio ap- 

i — I 

proach is going way beyond Ref. 9 : where several phc- 
nomenological parameters had to be introduced. Sec- 
ond, we study whether the Boltzmann equation approach 
with the LDA band dispersion as an input works well 
for this system. We will show that this approach gives 
results quantitatively similar to those of LDA+DMFT. 
Even though LiRli204 is a material very different from 
Na x Co02, having among others a much more compli- 
cated bandstructure, our analysis nontheless reveals that 
the origin of the large thermopower is similar: the "pud- 
ding mold" shape of the bands crossing the Fermi energy. 
This outcome was not prejudiced in our investigation and 
is quite surprising. We also discuss how electron doping 
could further increase the power factor of LiRh 2 04. 



II. METHOD 

As a first step we do a LDA calculation for LiRli2 04, 
using the linearized muffin tin orbital (LMTO) basis 
seti 1 ^, employing the experimental lattice constant a — 
8.46 and so-called x parameter x — 0.261 (which indi- 
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FIG. 1: Left panel: (Color online) Band dispersion of the 
effective 3-orbital Hamiltonian (solid line) and total LMTO 
band structure (dashed line) of L1RI12O4. Right panel: partial 
ai g and density of states for the model. 



cates the position of the oxygen sites) . 

From the LMTO bandstructure, we construct an effec- 
tive Hamiltonian (= H^ A ) by the projection onto Wan- 
nier functions^. Since the unit cell of L1RI12O4 contains 
four Rh atoms and each Rh atom has three ti g orbitals 
the size of the effective Hamiltonian is 12 x 12. A compar- 
ison of the band dispersion of this effective Hamiltonian 
with the total LDA band structure is shown in Fig. Q] 
In contrast to the case of Na^OoO^, not only the a% g - 
orbital but also the e^-orbitals have a substantial den- 
sity of states (see the right panel of Fig[TJ) at the Fermi 
level (Ep). Hence, we cannot extract a simpler effective 
Hamiltonian from the 12 x 12 Hamiltonian and need to 
keep all ti g orbitals in the following calculation. 

Next, we supplement the 3-orbital Hamiltonian by lo- 
cal intra- (U) and inter-orbital (U 1 ) Coulomb repulsions 
as well as by Hund's exchange (J; of Ising type), and 
solve it by DMFT 13 , using the quantum Monte Carlo 
(QMC) method. To get high-quality QMC data, we take 
~ 3.0 x 10 7 sweeps in the calculation. 

In the framework of DMFT, the Kubo formula for the 
Seebeck coefficient 
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where fc^ and e are the Boltzmann constant and unit 
charge, respectively, and 

/oo 
dwcjf x {w)f^)f{-u){^) n , (2) 
-oo 

= yY,Tr(v*(k)p(k,u;)v*Qi)p(k,u)). (3) 

k 

Here, p(k, u>) is the spectral function, i.e., the imag- 
inary part of the Green function G(k, w); ^ Q ^(k) = 



(kj3\(\/m)V x \ka) is the group velocity, f(u>) the Fermi- 
Dirac distribution function, and V the volume of the unit 
cell. 

As is carefully discussed in Ref. [I?], when the tight- 
binding basis is well localized in the real space, we 
can use the so-called Peirls approximation, u a/ g(k) = 
Vk-ff^ A (k) 14 . In this method, since we have an ana- 
lytical expression of H^ A (k) , the mesh for the momen- 
tum sum in Eq.Q can be arbitrarily dense. In most 
cases we took a 40x40x40 mesh, but in some cases also 
a 80x80x80 mesh for checking convergence. 

Usually, G(k, u) is calculated in DMFT(QMC) from 
the self energy which is obtained as a root from the 
local Green function G lmp (w), obtained in turn from the 
QMC data by the maximum entropy method (see, e.g., 
Refs. 18.19). However, this standard approach does not 
work well for the calculation of the Seebeck coefficient 
because of the following: Since 4> xx (w) only contributes 
to A m for \uj\ < ksT, we need for small \w\. For 

such frequencies E(w) is quite small (smaller than O.leV 
for \u\ < ksT, see below). As is pointed out in Reffl3. 
this smallness makes it difficult to calculate £(a>) reliably 
especially by a probability-based algorithms such as the 
maximum entropy method. 

Hence, in the present study, we calculate S(w) directly 
from T,(ioj), using both the Pade approximation and a 
polynomial fit. For the former, we apply the algorithm 
proposed in Ref. |2fj to the data with iuj £ [0, 45i]eV. For 
the latter, we fit S(?w) for iuj E [0, 4i]eV to 53«=0 c n(^ ,u> ) n 
by a standard least-squares fit. Since only the behavior 
at small |u;| is relevant for the Seebeck coefficient, we can 
expect the polynomial fit to give reasonable results. The 
Pade approximation might become problematic if poles 
are present in the vicinity of the real-w axis. However, 
as we will see below, the resulting S(w) for the present 
case does not show any anomalous behavior for small 
\uj\, which implies that the Pade approximation is not 
problematic. 

In Fig. El we plot £(w) for (77, V, J) = (3.1, 1.7, 0.7)eV 
which was estimated in Ref. [l8| and (3 = 1/fcgT =30, 34, 
40 eV _1 . For T ~ 300K, the main contribution stems 
from ui G [— 0.03, 0.03]eV. For these energies Pade ap- 
proximation and polynomial fit give similar results. Even 
though the agreement is not perfect, differences arc small, 
i.e., of 0(0.01) eV. Thus we employ E(w) of both, Pade 
approximation and polynomial fitting, in the following 
LDA+DMFT calculation of the Seebeck coefficient. The 
difference gives us an estimate for the accuracy of our 
calculation. 

Besides the LDA+DMFT study, we also performed cal- 
culations based on the Boltzmann equation. The Seebeck 
coefficient can be estimated by calculating 
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FIG. 2: (Color online) LDA+DMFT(QMC) self energy cal- 
culated by the Pade approximation (left) and a polynomial 
fit (right). 

Here, r is the relaxation time which we assume to be inde- 
pendent of k; e Q (k) are the eigenvalues of H^j/ L (k); and 

u a (k) the diagonal elements of Wv a ^{k)U , where U is 
the unitary transformation which diagonalizes H^p A (k). 
Note that K n can be roughly estimated as 

#o~5>d + «i). (6) 
#i~(teT)J>|-i£), (7) 

apart from a constant factor™. Here, is a summation 
over the states in the range of |e(fc)| < O(fcsT), and ua 
and ub are typical velocities for the states below and 
above the Fermi level, respectively. 

III. RESULTS 

In Fig. [3l we show the resulting Seebeck coefficient 
calculated by the LDA+DMFT method and the Boltz- 
mann equation approach. We also plot the result of the 
constant-r approximation for the Kubo formula, i.e., we 
assume E(w) = -1.0~ 3 i for Eqs. 

From Fig. [3l we see that (1) the Boltzmann equation 
and the constant-r approximation for the Kubo formula 



FIG. 3: (Color online) Thermopower calculated by the Boltz- 
mann equation approach and the constant-r method as well 
as by LDA+DMFT, using both the Pade approximation and 
a polynomial fit for the self energy. 



give almost the same result; (2) the constant-r approx- 
imation gives a larger thermopower than LDA+DMFT; 
(3) this smaller LDA+DMFT thermopower is closer to 
experiment^, already for (U, U', J) = (3.1,1.7,0.7)eVi^ but 
even more so for somewhat smaller values of the Coulomb 
interaction. 

Point (1) demonstrates that the calculation via Eq. ((T|) 
is working well, if £(w) is correct. Point (2) can be un- 
derstood from the behavior of — Imll(u;): Fig. [2l shows 
that -ImE(w) calculated by the LDA+DMFT method 
is large for negative oj but small for positive 10 (indepen- 
dently of Pade approximation and polynomial fit). This 
means that, in contrast to the constant-r approximation, 
the actual life time of quasi-holes is longer than that for 
quasi-particles. Therefore, the contribution of the quasi- 
holes (particles) to 4> xx in Eq. ([!} becomes larger (smaller) 
in the LDA+DMFT calculation, and consequently the 
first moment, A\, becomes smaller. Here, it should be 
noted that the constant-r approximation does not cor- 
respond to the limit of U = U' = J = 0, since this 
asymmetry of life time exists even in the weak coupling 
limit. This is the reason why the results of LDA+DMFT 
move away from those of the constant-r approximation 
as U, U' and J are decreased. 
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FIG. 4: (Color online) Group velocity squared (u k ) 
along different directions of the fisrt Brillouin zone for 
Rh+ 3 ' 5 (LiRh 2 4 ; upper panel) and Rh +3 ' 08 (electron-doped 
LiRli204; lower panel), k point above the Fermi energy Ef 
are shown in yellow, those below Ef in red. 
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FIG. 5: (Color online) Power factor (normalized by its value 
at Rh valence=+3.5 ) and thermopower (inset) as a function 
of the valence of Rh, calculated by the Boltzmann equation. 



As for point (3), we would like to note that the corre- 
lations renormalize the bandwidth. This renormalization 
is calculated microscopically here whereas it has been ad- 
justed to the angle-resolved photoemission spectrum in 
Ref. & 



IV. DISCUSSION 

While there are some differences between Boltzmann 
equation approach and LDA+DMFT, the results are still 
very similar, even quantitatively. Hence, we may expect 
that, in the present case, the Boltzmann equation can be 
used as a convenient tool to analyze the mechanism of 
the large thermopower, or even to design more efficient 
thermoelectric materials. 

Let us first examine whether the mechanism proposed 
for Na x Co02 in Ref. |9| can work also in LiRli2 04. If the 
valence band has a peculiar shape of dispersion which is 
dispersive below the Ep but somewhat flat above (the 
so-called "pudding-mold" type), K\ in Eq. ([7]) becomes 
large, since the group velocity above Ep {u 2 A ) is much 
larger than the one below Ep {u 2 B ) in this case. This 
is the basic idea of Ref. how to realize a large ther- 
mopower and a low resistivity at the same time^i. In the 
top panel of Fig. QJ we plot the group velocity squared for 
LiRh204 within the energy window of |e — Ep\ < 3k B T 
at T ~ 300K. We see that u\ is indeed larger than u 2 B , 
confirming this mechanism. We note here that although 
the Rh valence is +3.5 in LiRli2 04, the degeneracy of d xy , 
d yz and d zx orbitals in the cubic phase makes the number 
of holes per band small, resulting in a situation similar to 
Na x Co02 with the Co valence smaller than +3.5. This 
view is consistent with the experimental fact that the 
thermopower is suppressed in the tetragonal phase be- 
low 230K, where the degeneracy is lifted. 5 

However, we also see that the (squared) group velocity 



above Ep is still large for some fc-points. In fact, for 
LiRh204, there are two pudding mold bands. For the 
Rh valence of +3.5, Ep lies near the bending point of 
one of the pudding mold bands, but also cuts through 
the dispersive portion of the other (see the upper inset 
of Fig. [4]). The former enhances the thermopower, while 
the latter suppresses it. This might be the reason why 
LiRh2C>4 is not such a good thermoelectric material as 
Na a; Co02. 

To enhance the thermopower in LiRh204, we suggest 
the following possibility: If we dope electrons to this sys- 
tem, one of the pudding mold bands will be brought com- 
pletely below Ep (see the lower inset of Fig. [4|. The 
second panel of Fig. |4] shows the group velocity squared 
for such a doped system with Rh valence +3.08. In this 
case, the group velocity squared above Ep is small for 
the entire Brillouin zone. 

To confirm this idea, we calculate the thermopower 
and the power factor (normalized by its value at Rh 
valcncc=+3.5) for various Rh valences by means of the 
Boltzmann equation approach^, see Fig. [5] The results 
indicate a maximal power factor (= S 2 / p oc Kf/Kg) at 
a valency of +3.08, where it is almost four times larger 
than for L1RI12O4. While the orbital degeneracy of d xy , 
dy Z and d zx already plays a crucial role to make Ep be 
higher than those of single-orbital systems- realizing this 
situation experimentally is an interesting challenge which 
seems to be feasible. 



V. ACKNOWLEDGMENT 

We would like to thank H. Takagi and Y. Okamoto 
for fruitful discussions and providing the experimental 
data shown in FigJ3] Numerical calculations were per- 
formed at the facilities of the Supercomputer center, 



5 



ISSP, University of Tokyo. This work was supported grant 19019012,19014022, 19051016 and Russian Foun- 
by Grants-in-Aid for Scientific Research (MEXT Japan) dation for Basic Research (RFBR) grant 07-02-00041. 



Present address: Department of applied physics, Univer- 
sity of Tokyo, Tokyo 113-8656, Japan. 

1 G. D. Mahan, Solid State Physics 51, 81 (1997). G. D. 
Mahan et al, Physics Today, March 1997, p.42. 

2 A. Bentien, et al. Europhys. Lett. 80, 39901 (2007). 

3 S. Paschen, Thermoelectric aspects of strongly correlated 
electron systems, in CRC Handbook of Thermoelectrics, 
Ch. 15, (ed. D. M. Rowe, CRC Press), Boca Raton, 2005; 
A. Bentien et al, Europhys. Lett. 80, 17008 (2007). 

4 I. Terasaki, Y. Sasago and K. Uchinokura, Phys. Rev. B 
56, R12685(1997) 

5 Y. Okamoto et al, Phys. Rev. Lett. 101, 086404(2008). 

6 W. Koshibae et al, Phys. Rev. B 62 6869 (2000). 

7 D.J. Singh: Phys. Rev. B 61 (2000) 13397. 

8 GB. Wilson-Short et al: Phys. Rev. B 75 (2007) 035121. 

9 K. Kuroki and R. Arita, J. Phys. Soc. Jpn., 76, 083707 
(2007). 

I. Terasaki, JPSJ Online-News and Comments [Oct. 10, 
2007]. 

1 Y. Ishida et al., J. Phys. Soc. Jpn. 76, 103709 (2007). 

2 V. I. Anisimov et al, J. Phys. Condens. Matter 9, 7359 
(1997), A. I. Lichtenstein and M. I. Katsnelson, Phys. Rev. 
B 57, 6884 (1998); K. Held et al, phys. stat. sol. (b) 243, 
2599 (2006); G. Kotliar et al, Rev. Mod. Phys. 78, 865 
(2006); K. Held, Adv. Phys. 56 829 (2007). 

3 W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 
(1989); A. Georges and G. Kotliar, Phys. Rev. B, 45, 6479 
(1992); A. Georges et al, Rev. Mod. Phys. 68, 13 (1996). 

4 V. S. Oudovenko et al, Phys. Rev. B 73, 035120 (2006). 

5 O. K. Andersen, Phys. Rev. B 12, 3060 (1975); O. Gun- 



narsson, O. Jepsen, and O. K. Andersen, Phys. Rev. B 27, 
7144 (1983) 

16 V. I. Anisimov et al, Phys. Rev. B 71, 125119 (2005). 

17 I. Paul and G. Kotliar, Phys. Rev. B 67, 115131 (2003). 

18 Z. Pchelkina et al, Phys. Rev. B 75, 035122 (2007). 

19 I. Nekrasov et al, Phys. Rev. B 73, 155112 (2006). 

20 H. J. Vidberg and J. W. Serene, J. Low Temp. Phys. 29, 
179 (1977) 

21 Note that this idea is totally different from the conven- 
tional idea to realize large thermopower exploiting an 
asymmetry of density of states (DOS) around the Fermi 
level. While DOS roughly scales as (l/u(k)), the quantity 
which plays crucial role here is (u 2 (k)). These quantities 
behave quite differently in general. For example, for the 
case of Na2;Co02, if we calculate the former taking into 
account the dispersion along the z-axis, we see that it has 
no special structure around the Fermi level. However, on 
the other hand, the latter has a strong asymmetry, which 
we believe to be the origin of large thermopower in this 
material. 

22 Note that this calculation neglects the energy and filling 
dependence of r, which should be present and affect p and 
the power factor. Nevertheless, since the effect of the fre- 
quency dependence of r is not so serious in the calcula- 
tion of the Seebeck coefficient (Fig. [3]), we believe that the 
present estimation of the power factor at least gives the 
correct tendency of the filling dependence, especially if the 
Rh valence does not change drastically. 



